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^ ' Abstract 



Reaction equations of homogeneously mixed pollutants in the atmosphere can lead 
to non-stationary periodic solutions. It is important to know in which respect these 
[^ ■ solutions are modified under the influence of the atmospheric currents. We inves- 

^P ! tigate this question in a very simple model: The reaction equations are modeled 

by the equations of the brusselator and the currents are represented by an isolated 
vortex. In the limit of high vortex currents we find again the homogeneous solutions 
t^^ , whereas for smaller currents complicated spatial and temporal patterns emerge. The 

■»^ ' role of the diffusion as a singular perturbation is investigated. 

C/3 ' 

O 

• 1— I 
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'q,: 1 INTRODUCTION 



The concentrations of pollutants in the atmosphere depend on the atmospheric 
5^ I currents and the reaction equations between the pollutants. (Moreover they 

depend on cloud formation, humidity, ice etc, influences that will not be dis- 
cussed here.) 

When dealing with this problem the currents are usually replaced by their 
average and mixing is modeled by introducing 'turbulent diffusion' [1] [2] [3] [4] . 
The problem of such approximations is that the chemical reactions depend 
on the local concentrations and not on averaged ones. A further disadvantage 
of this scheme is that the inserted turbulent diffusion is orders of magnitudes 
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higher than the molecular diffusion. This is questionable because diffusion is 
a singular [5] [6] and sensitive perturbation as will be shown below. 

If the atmospheric currents lead instantaneously to a homogeneous mixing 
of the pollutants, all concentrations are obtained from the chemical reaction 
equations alone. The solutions of these reaction equations need not approach 
a time-independent limit (fixed point for the dynamical system) but periodic 
fluctuations are possible as has recently been shown in a model system contain- 
ing six pollutants and two pollutant sources[7][8][9]. A particular feature was 
the result that pollutant concentrations may change by an order of magnitude 
within a few days. 

The solutions of the reaction equations become perturbed and then modified 
as soon as the atmospheric currents do not lead to a complete mixing any 
more. This situation could be modeled by adding to the reaction equations 
the above mentioned averaged currents plus a turbulent diffusion term that is 
sufficiently high. But according to the remarks made earlier a more appropri- 
ate approach is to assume that the diffusion is weak and to incorporate the 
transport equations in the reaction equations directly. 

The temporal evolution of concentrations of n reactants are described by the 
react ion- advect ion- diffusion equat ions : 

d 

-^Ci = fi{ci..Cn,yi,t) -Y{:si,t) -Vci + eiAci, i = l,..,n, (1) 

where the functions /« describe the chemical reactions and pollutant sources. 
These functions depend explicitly on the spatial coordinate x if the reactant 
i is produced by a localized source. The second term represents the advection 
of pollutants by the velocity field v(x, t) and the last one is the diffusion 
term. The coefficient e is small and in a first approximation we neglect the 
diffusion. Later on we will investigate the addition of diffusion in our model 
very carefully. 

In general the equations eq.(l) are rather involved. In this paper we choose the 
simplest possible system describing nontrivial reactions between constituents 
subjected to a simple current which is a solution of Euler's equation[10]. The 
reaction equations are taken from the brusselator[ll]. Depending on the pa- 
rameters the concentrations converge to a fixed point or limit cycle in the 
homogeneous case. Typically periodic solutions with sharp peaks occur. They 
are analogous to those obtained in the model of [7] where the chemistry of 
tracer constituents in the troposphere were described. 

For maintaining the chemical reactions a 'pollutant source' is required. We 
choose a point source pouring out pollutants of sort X placed in the veloc- 
ity field of a two-dimensional isolated point vortex. The pollutant X decays, 



partly into harmless substance E partly into a second constituent Y that au- 
tocatalytically reacts with X again according to the well known reactions of 
the brusselator 



X^Y (2) 

2X + Y^3X 



Since molecular diffusion does not lead to a spreading on a macroscopic scale 
the vortex distributes the pollutants along the circular streamline containing 
the point source and thus the model becomes one dimensional. We note that 
this advection-reaction problem can be seen as an extension of the simple 
brusselator containing it as a limit limiting case for high vortex strengths (or 
when the source is close to the vortex center). 

Next we discuss the distribution of concentrations on the circle. When the 
time T, needed for one circulation is small, the solutions become very similar 
to the reaction equations with a homogeneously distributed source. On the 
other hand, if T is larger than a threshold value for an observer moving with 
the fluid, the distribution becomes periodic with period T, whereas for an 
observer at rest the distribution is stationary. When T decreases the distri- 
butions become rather different: Time dependent solutions with period nT, 
{n = 2, 3...), quasiperiodic and chaotic solutions are detected. All these lead to 
time dependent distributions that are infinitely degenerate on the circle and 
therefore depend on the initial distribution. Even if that is smooth, steps in 
the distributions occur, in the chaotic case on each scale. Therefore diffusion 
is a singular perturbation and switching it on along the streamline leads to 
drastic effects. For the period 2T case, the addition of diffusion can be un- 
derstood completely. We find two time scales. The first, which is of the order 
of iQ leads to coarse graining on the order of y/e, the second which is of the 
order of e^^^, P = 0{1) removes all the degeneracies of the solutions leading 
to a distribution that is smooth apart from one step. This step is intrinsic and 
would appear for a source of finite size along the streamline as well. 

For parameters, that without diffusion lead to periodic solutions of higher 
period we observe pattern formation and spatiotemporal chaos. All these so- 
lutions have nothing in common with the case of homogeneous mixing we 
started with. 

In section II we present the model, in section III we discuss its properties 
without diffusion, in section IV we concentrate on the role of diffusion. The 
conclusion ends the paper. 



In this paper we use non-dimensional units. 



2 THE MODEL 



The reaction processes, eq.(2), lead to the well known reaction equations, 
which in non-dimensional form read as 



ci = c?C2-(l + 6)ci (3) 

C2 = bCi - clc2 

Here ci and C2 are the concentrations of constituent X and Y respectively. 
The parameter b presents the ratio between the decay rate of X into harmless 
substances and the decay rate of X into pollutant Y. 

Without a source term the concentrations tend to zero. In our model we 
assume for simplicity a point like pollutant source that is fixed, but encounters 
the fluid field of a two-dimensional isolated vortex. This vortex produces a 
circular flow with velocity of modulus T/r at distance r from the center, F 
being the strength of the vortex. This problem is one-dimensional since the 
chemical reactions can only take place along the circular streamline containing 
the source and can be mapped to the unit interval with periodic boundary 
conditions. Thus the flow is completely represented by one parameter, the 
dimensionless velocity v. We obtain the following combined reaction transport 
equations: 



d d d"^ 

— ci = s6{x) + clc2 - (1 + b)ci - v—ci + e^ci (4) 

— C2 = bci - qc2 - V--C2 + e—^C2 
at ox ox^ 

where s is the strength of the source located at a; = 0. We allow here for a 
very small amount of diffusion with diffusion constant e. 

Apart from the diffusion constant e the equations contain three relevant pa- 
rameters, the decay ratio 6, the source strength s and the velocity v. 

Eq.(4) represent the Eulerian description of an observer at rest. If we change 
to the Lagrangian description (i.e. to an observer moving with the fluid) eq.(4) 
is transformed into 



—ci = s ^ 5(x + t;t + n) + clc2 - (1 + b)ci + e^r^ci (5) 

ot „^g ox 



-C2 = 6c,-c,C2 + e- 
x = x — vt mod 1 



The solutions to this model will be discussed in the following two sections. 



3 PROPERTIES OF THE MODEL WITHOUT DIFFUSION 



First let us discuss the transport reaction equations, eq.(4) and eq.(5) with 
diffusion switched off: 

For e = the coordinate x can be interpreted as a parameter equivalent to a 
time translation x/v in the driving term. Therefore in the moving frame we 
have to solve for any x the ordinary differential equations 



d °° 

— ci = s ^ 5{x + vt + n) + c\c2 - (1 + h)ci (6) 

^* n=0 

d , 2 

--Co = OCi — Ci Co 

dt ^ 

The spatial distribution of the concentrations is completely determined by the 
solution of the above equation and the initial distribution. Obviously we have 
to take into consideration the time shift of the driving for different points. 

In eq.(6) the difference to the 'usual' brusselator is the periodic 5 function 
time dependence of the source. Thus we have a periodically driven (kicked) 
brusselator with one extra parameter, the period of the driving T = 1/v. 
The periodically driven brusselator has been investigated in different contexts 
considering a constant plus a sinusoidal or delta function time dependence of 
the source [13] [14] [15]. In the T ^ limit the normal brusselator is recovered, 
i. e. very frequent injections correspond to an almost uniform source. In this 
limit the parameter plane s — h can be divided into two regions (fig.l): For 
values corresponding to higher source strengths the concentrations converge 
to the fixed point c\ = s, C2 = h/ s. As s is decreased, the fixed point becomes 
unstable and a Hopf bifurcation occurs along the curve s = \Jh + 1 forming 
the boundary between the two regions. Below this curve the system converges 
to a limit cycle, i. e. the concentrations oscillate periodically, (fig.2.) 

As the periodic driving is switched on < T <^ 1, a periodic pulsation with 
period T of the concentrations appears. Moreover the initially two dimen- 
sional phase space becomes three dimensional by including the cyclic variable 
tjT mod 1 due to the driving. Thus the dimensionality of the attractor in- 
creases as well, and the original fixed point turns to a limit cycle, representing 
a periodic time dependence, but still remains a fixed point of the stroboscopic 
map. Similarly, the original limit cycle becomes either a torus corresponding 
to a quasiperiodic dynamics with one of the periods equal to T (fig.3) or a 
periodic orbit with large period. This is expected from the characteristic fea- 
tures of periodically driven oscillators. There exist resonant regions for driving 



frequencies close to their natural frequency multiplied by a rational number. 
These resonant regions appear here below the Hopf bifurcation curve. In the 
parameter space of these regions Arnold tongues are detected analogous to 
those of the so called 'circle map' [16] (fig.4,5). As another effect of the driv- 
ing the Hopf bifurcation curve moves to smaller values of s as T increases 
(fig.l). Since the dynamics is given by a set of two non-autonomous ordinary 
differential equations chaotic behavior is also possible for some values of the 
parameters leading to a strange attractor in the stroboscopic section (fig. 6). 

Next we consider the distribution of the concentrations due to the properties 
of the driven system. Let us discuss the typical cases: 
a)periodic time dependence with period T: 

This behavior occurs for large T, i.e. as the velocity of the flow becomes small. 
In this case the concentrations oscillate and the phase of the oscillations is 
determined by the phase of the driving t/T mod 1. Thus the final state does 
not depend on the initial conditions. In the moving frame the only difference 
in the periodic time dependence at different points of the flow is a time lag 
xT. 

c{x,t) =c{0,t + xT). (7) 

If we move back to the standing frame (a; — * a; — vt) we obtain asymptotically 
a stationary state c{x, t) = c(0, xT) i.e. the concentrations converge to a x- 
dependent function constant in time. 

f3) periodic oscillations with period uT: 

This behavior corresponds to the resonant regions. In this case the concentra- 
tions can take n different values for a given phase of the driving depending on 
the initial conditions. 

c{x,t) = {c{0,t + iT + xT)}, i = 0...n-l, (8) 



The boundary between the basins of attraction of the n branches of the solu- 
tion is a twisted (Mobius-like) surface so that the basin of attraction of branch 
i becomes the basin of attraction of branch i + 1 mod n after one period T. 
Thus any smooth initial condition must have at least one intersection with this 
surface. At this point the concentrations converge to two different branches so 
a discontinuity appears in the spatial distribution of the concentrations (fig. 7). 
Note that this step is not a consequence of the delta function in eq.(6) but is 
due to geometrical constraints. An initially random distribution can lead to a 
completely staggered distribution whose envelopes are the n branches of the 
solution. In the moving frame the steps remain at the same position x. 

7) quasiperiodic time- dependence: 

This is present in a region below the Hopf bifurcation curve between the 



resonances and is pronounced for small T because with increasing T the size 
of the region below the Hopf curve shrinks and at the same time the resonant 
islands grow in size. This case corresponds to a motion on a torus in the phase 
space. The dynamics can be characterized by two cyclic angle-like variables, 
one of them is the phase of the driving and the other one depends smoothly on 
the initial concentrations. c{x,t) = c{x,t + T{x)) = c{t + T{x) +xT) Therefore 
an initially smooth distribution remains smooth in x for all times (except at 
the initial position of the source where the time lag of the driving by T leads 
to a discontinuity.) 

5) chaotic time- dependence: 

In this case the time dependence is very sensitive to the initial conditions and 
thus the distribution becomes irregular on each scale regardless how smooth 
the initial distribution may have been (fig.9a). 



4 THE ROLE OF DIFFUSION 



Without diffusion the final distributions (except those with the period T) have 
infinite degeneracy due to an arbitrary uneven[3] number of steps. Therefore 
diffusion is a singular perturbation that has significant consequences for the 
system as small as e may be. For the following computations we used the 
Crank-Nicholson scheme combined with operator splitting [12]. 

We discuss here the simplest nontrivial case of eq.(8) first, which occurs for 
period 2T. We denote with a — h (-1 — ) step an 'upward' ('downward') steep 
increase (decrease) of the concentration, but exclude the strong increase of 
Ci at the location of the source. If the diffusion is small enough we can treat 
a step as isolated (for a very long time). Due to diffusion the step will move 
with a drift velocity, cf (fig. 7). Scaling and symmetry arguments suggest that it 
should be proportional to a higher power of ^/e and in fact we find numerically 
a dependence oc e. The important point, however, is that, averaged over 2T, 
each isolated step moves with the same drift velocity. In fact after time T 
a — h step becomes a H — step and vice versa. What we expect then as the 
essential ingredient of eq.(8) is that / tries to enforce solutions with period 
of 2T. In appendix A we have derived a simple function f that has just this 
property anc? makes it possible to treat eq.(8) analytically. Then we find: first, 
for times 

Ti = 0{1) (9) 



^ we do not count the strong increase of the ci concentration due to the 5 function 
shape as a step. 



all steps with distance of 0{^/e) vanish. During this time the diffusion does 
nothing but a coarse graining. Second, over a period of about 

T2 = e^/v^, (3 = 0{l) (10) 



all other steps are effected. The diffusion removes the degeneracies by and by, 
until a final state emerges that has no steps at all, besides the generic one that 
cannot be removed. This state has global stability in our model. Numerically 
we find the same phenomena for f of the brusselator, cf. fig. 7b. 

The effects described here occur quite independently from how small e is again 
demonstrating that the diffusion is a singular perturbation. On the other hand 
T2 depends exponentially on 1/^/e. When a further perturbation has to be 
added acting on a time scale r we expect quite different situations depending 
on whether T2 > r or T2 < r. This means that the effect of such a pertur- 
bation depends sensitively on ^/e which shows that introducing diffusion as a 
parameter - as has been done by introducing 'turbulent diffusion' - is quite a 
dangerous approximation. 

We expect even more complicated properties of the concentrations having 
higher periods (in absence of diffusion). There are two reasons for that: i) if the 
period is nT the system has at any location n— 1 choices for the height of a step, 
ii) the steps are no longer equivalent but are separated in classes and only steps 
within the same class change into each other and therefore move with the same 
mean drift velocity Vd- Indeed, the effect of diffusion on the periodic solution 
can be very significant in some cases by leading to a complicated irregular 
behavior of the system in space and time. As can be seen from fig.8, inside 
the chaotic concentration field coherent regions with regular periodic time 
dependence appear and disappear continuously. This kind of spatiotemporal 
intermittency has been observed in different extended systems[17])[18]) (e.g. 
in case of coupled maps[19]) [20] [21]. If one starts the simulation with a smooth 
initial distribution first at least the intrinsic step appears as described above. 
The pertubation of the periodic solution around the discontinuity leads to a 
chaotic time dependence which due to the diffusive coupling spreads over the 
whole system. Such behavior can be observed for parameters which lie in the 
vicinity of the chaotic regimes in the e = case. The solution appears already 
for very small e demonstrating again that diffusion is a singular perturbation. 

In case of quasiperiodic local behavior, instead of a finite number of discrete 
branches, a continuous set of solutions exists filling the torus in the phase 
space. Thus the discontinuity present in the case without diffusion is easily 
removed by an arbitrarily weak diffusion leading to coherent quasiperiodic 
oscillations of the whole system. 

When the parameters correspond to chaotic local dynamics, diffusion tends to 



form correlated regions of finite extent in space and time (fig- 9) similar to the 
case above. As e is increased, the local dynamics becomes completely regular 
with a frozen irregular distribution in space which certainly depends on the 
initial distribution. 



5 Conclusion 



High peaks can appear in periodic solutions of chemical reaction equations 
in which the constituents are homogeneously mixed tracer gases of the atmo- 
sphere. However, depending on the motion of the fiuid the mixing need not 
be homogeneous at all, and the question arises how these solutions will then 
change. 

In this paper we investigate this question for a simple model, the brusselator 
with pointlike source in a one-vortex fiow. Simple as the model appears, it 
already demonstrates the strong modifications occurring as soon as we move 
away from the homogeneous situation. One observes this when computing the 
concentration distribution along the (closed) streamline in which the source 
is located. As a function of time we detect solutions that are very similar 
to those of the homogeneous case. This happens as long as the period T of 
the fiow is small. Furthermore we find solutions with period nT, moreover 
quasiperiodic and chaotic ones. All these solutions, except that with period 
T are infinitely degenerate and therefore depend on the initial distribution. 
Even if that is smooth, the distributions can have asymptotically an arbitrary 
(uneven) number of discontinuities, in the chaotic case on each scale. 

In such situations diffusion is a singular perturbation and switching on arbi- 
trary small diffusion along the streamline has two effects: First after a time of 
C(l) it leads to a 'coarse graining' of the distribution on a space scale oc ^/e 
where e is the strength of the diffusion. Second on a time scale oc e^^^, 
a = 0{1) it removes all discontinuities but one for solutions which have (with- 
out diffusion) period 2T. This shows that the solutions depend sensitively on 
^/e. For parameters that lead (without diffusion) to solutions of higher period, 
quasiperiodic solutions and chaotic ones were observed, moreover pattern for- 
mation and spatio temporal chaos. All those solutions have nothing in common 
with the case of homogeneous mixing we started with. 

Although this one-dimensional model is far from being a realistic represen- 
tation of the chemistry and transport in the atmosphere, it shows that even 
a trivial non-turbulent flow interacting with a simple regular chemical dy- 
namics of just two reactants can lead to a complex irregular behavior of the 
concentration fields. 
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A appendix 

In this appendix we derive the properties of eq.(l) for our model assuming 
that without diffusion the solution has a 2T period in the moving system. We 
have 

9iC = f (c, X, t) + e(9?c (A.l) 

with the periodic boundary conditions 

c{x,t) = c{x+l,t) (A.2) 

dsic{x,t) = dsc{x+ l,t) 

This equation holds true in the frame moving with a velocity v = 1/T. 

The use of c can be awkward since the components have to be positive. There- 
fore we write 

n = c + const (A. 3) 

and get the equation for n 

dtn = g(n, X, t) + edln 

with (A.4) 

g(n, x, t) = f (const + n, x, t) 

and the periodic boundary conditions 

n{x,t)=n{x + l,t) (A. 5) 

ds:n{x,t)=dsn{x + l,t) 

Without diffusion n moves exponentially fast to its asymptotic limit n*^°^ hav- 
ing the properties 

n(°)(x,t) = n(°)(a:,t + 2T) 

and either 
n(°)(x,t) = n(°)(t-xT) (A.6) 

or 
n(°)(a;,t) = n^") (T + t - xT) 
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n*^°^ and its properties remain important also if diffusion is switched on since 
g can be expanded around n^^\ 

To understand the physics of eq.(A.4) with conditions eq.(A.6) we construct 
a simple model for the function g in three steps: 



(0) 



step 1: we introduce a very simple n 



nf\x,t) = ^{ae'^^'/^-''^} (A.7) 

and we represent the two dimensional vectors by complex numbers, 
step 2: we use the ansatz 

n(a;,t) = e'"(*/^-^)-m(x,t) (A.8) 

and get the pde for m 

dtm = g — medxTn — Ti^em + ed^m 

with (A.9) 

The terms medxfn and vr^em are of higher order in ^/e and will be left out for 
simplicity. The boundary conditions of eq.(A.5) are replaced by 

m(x,t) = -m(x + l,t) (A.IO) 

dxTn{x, t) = —dxTn{x + 1, t) 

step 3: we construct a simple g. Because of eq.(A.7) and eq.(A.8) m*^°^ can 
take two values only, 

m(°) = ±a (A. 11) 

and a can be chosen to be real and positive. When m is in the neighborhood 
of m^^^ g can be expanded and we obtain 

g = —a{x, t){m — a) + ... 

or (A.12) 

g = —a{x, t){rn + a) + ... 
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For a{x, t) we insert a real positive constant^. The linear approximation of g 
is of course incorrect if m is not close to ±a. A nonlinear ity is added simply 
by the prescription 

—aim — a) for \ m — a \<\ m + a 
9 [ 

-a{m + a) else 



Thus we get the partial differential equation: 

dtm = g + edlm 

with (A. 13) 

—a{m — a) for ^m > 

— Q;(m + a) for ^m < 

Boundary conditions are given by eq.(A.lO). The connection between n and 
m is given by eq.(A.8) and real and imaginary part of n are the components 
of n. 



properties of the solutions of eq.(A.13) 



I) Diffusion switched off, i.e. e = 0, 

m consists asymptotically of an uneven number of steps with values ±a. The 
number of steps can be arbitrarily high and is determined exclusively by the 
initial distribution of m. 



II) Diffusion switched on, i.e. e > 0, 

1) Isolated step 

we assume that there is a constant velocity ^/ew with which the step is moving. 
Transforming to new coordinates y with 

x = y + y/ewt (A. 14) 



a could be a complex constant as well as long as the real part is positive. 
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and assuming that the step occurs at y = the two equations are to be solved: 

= —a{m + a) + y/ewdyin + edyin, y < 
= —a{m — a) + ^/ewdyin + edyiJi, y > 

because of the boundary conditions for the isolated step 

m(— oo) = —a 
m(oo) = a 

the solution is 

m_(i/) = Ae''^/^-a, 7 = {l/2){-w + V4:a + w^), y<0 (A.15) 

m+{y) = Be^y^^' + a, 7 = (l/2)(-w - V4:a + w^), y>0 

with the boundary condition 
m_(0) =m+(0) 

Because of eq.(A.13) there is the further condition 

sR{m_(0)} = 

Therefore A = a and B = —a, m is real and the condition for w is obtained 
from 

07 = —07 
which means 

w = 

2) two interacting steps isolated from the rest 



Let the — h step be left, the H — step be right. Both steps move because of 
interacting with each other and we assume that the interaction changes speed 
and shape of the steps only slowly (the distance 2xd between them decreases 
of course). 

First we rescale to avoid the e dependence 

?^^ (A.16) 
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Next we transform into a coordinate system moving with the — h step, whose 
position is at 0. we obtain 

^ = V + Wdt (A. 17) 



and 



= —a{m^ + a) + Wddr^m^ + d^rri-, rj < (A. 18) 

= —a{m^ — a) + Wddr^m^ + i9^m+, r] > 

Boundary conditions: 

m_(— oo) = —a 

m_(0) = m+(0) (A.19) 

m'_{0)=m'^{0) 

Furthermore the presence of the H — step is taken care of by the condition 

m'4^d)=0 (A.20) 

and we have the constraint (cf eq.(A.13)) 

gfj{m_(0)} = (A.21) 

Then we get with an exponential ansatz (cf eq.(A.15) 

m_ {r]) = Ae''' -a, r]<0, (A.22) 

m+{r]) = Be""' + Ce^" + a, < r] < ^a 

Again m can be chosen to be real and the conditions eq.(A.19), eq.(A.20) and 
eq.(A.21) yield 

A- a = B + C + a 

^7 = ^7 + C7 (A.23) 

= B^e''^" + C-fe^^^ 
A-a = 

From these equations we get Wd (neglecting all terms w'^ and higher) 

Wd ~ 2y/^e~^^^^ (A.24) 
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which is correct for 

V^U > 1 (A.25) 

From scahng arguments we infer that eq.(A.25) gives the correct order of 
magnitude for ^/aC,d < 1- 

One can use eq.(A.24) also to prove that there is no stationary state. If it were 
all the equations were exact and in particular eq.(A.24) which in turn would 
be a contradiction. 

Now we compute the lifetime ti of a step which is in the original coordinates 

t, = -^(e2v^^''/v^-l) (A.26) 



3) n interacting steps 



To treat this problem we take into account the interaction between nearest 
neighbors only ( the interaction between next nearest neighbors is exponen- 
tially small compared to the interaction between the nearest neighbors). Then 
it is sufficient to look into the problem of one step between two other steps. 
We approximate the interaction again by boundary conditions and obtain two 
conditions of the form eq.(A.20). Doing an analogous computation with the 
same approximations we obtain for the velocity of the step 

Wd ^ 2v^(e-2V^^<^+/^ - e-2v^^d-/v^) (A.27) 



Here 2xd+ {^Xd^) is the distance to the right (left) step. From this result we 
conclude that all states with more than one step will be unstable since two 
neighboring steps will annihilate each other. 
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Fig. 1. Curves corresponding to the Hopf bifurcation in the parameter plane s — b 
for different values of the period T. 

Fig. 2. Constant in time and periodic behavior of the concentrations ci and C2 for 
the unforced brusselator (T = 0). The parameters are s = 2.5, b = 3.0 and s = 1.0, 
b = 3.0, respectively. 

Fig. 3. Quasiperiodic time dependence of the concentrations ci and C2 (a), and 
stroboscopic section (b) for s = 1.0, b = 3.0 and T = 1.0. 

Fig. 4. Stroboscopic plot of C2 in function of T for s = 1.9 and b = 7.7. The 
Hopf bifurcation occurs around T = 1.03 and there are resonant windows inside the 
quasiperiodic region labelled by the ratio of the two periods. 

Fig. 5. Periodic(blank) and quasiperiodic(gray) regions in a section of the parameter 
space for T = 1.0. The behavior of the system was identified by calculating the 
leading Lyapunov exponent which is smaller then —0.0025 for the blank region. 

Fig. 6. Chaotic time dependence of the concentrations ci and C2 and the strobo- 
scopic section of the strange attractor. The parameters are s = 1.2, b = 7.0 and 
T = 1.36. 

Fig. 7. Spatiotemporal plot of the concentrations ci (left) and C2 (right) along the 
streamlines in the co-moving frame represented on a grayscale, so that concentra- 
tions increase from black to white. The simulation was started with both concentra- 
tions equal to zero and the initial position of the source is at a; = 0.2. Parameters 
are s= 1.0 6= 5.0 and T = 1.7 In case (a) e = and a non-moving discontinuity 
is present at x = 0.2. When diffusion is switched on e = 0.001 the discontinuity 
becomes rounded and moves (to the left in this case) along the streamline. 

Fig. 8. Stroboscopic spatiotemporal plot of concentration C2 for parameters s = 0.8, 
b = 6.0 and T = 1.85, that correspond to a periodic behavior with period 3T when 
diffusion is neglected. Here e = 2 • 10^^ that leads to an irregular spatiotemporal 
dynamics. 

Fig. 9. Stroboscopic spatiotemporal plots of concentration C2 for s = 0.8, b = 6.0 and 
T = 1.89. This parameters correspond to a chaotic local dynamics when diffusion is 
not considered. We assumed that the initial concentrations are randomly distributed 
in a small interval [0, 0.0001] for both constituents. The diffusion coefficient is e = 
(a), e = 1.5 • 10~^ (b) and e = 2 • 10"^ (c), respectively. 
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